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Abstract 

In this work we formulate the reaction theory of parity violation in compound 
nuclear states using Feshbach's projection operator formalism. We derive in 

> 

this framework a complete set of terms that contribute to the longitudinal 

o 

O 

also discuss the parity violating spreading width resulting from this formal- 
ism. We then use the above formalism to derive expressions which hold in the 
. case when the doorway state approximation is introduced. In applying the 

£3 ■ theory we limit ourselves in this work to the case when the parity violating 



X 



asymmetry measured in experiments with polarized epithermal neutrons. We 



(PV) potential and the strong interaction are one-body. In this approxima- 
tion, using as the doorway the giant spin-dipole resonance and employing well 
known optical potentials and a time-reversal even, parity odd one-body in- 
teraction we calculate or estimate the terms we derived. Among others, we 
evaluate numerically a new "direct" term and the channel coupling term first 
derived by Lewenkopf and Weidenmiiller, this in addition to the evaluation of 
the main "compound" term. In our calculations we explicitly orthogonalize 
the continuum and bound wave functions. We find the effects of orthogonal- 
ization to be very important. Our conclusion is that the present one-body 
theory cannot explain the average longitudinal asymmetry found in the recent 
polarized neutron experiments. We also confirm the discrepancy, first pointed 



out by Auerbach and Bowman, that emerges, between the calculated average 

asymmetry and the parity violating spreading width, when distant doorways 

are used in the theory. 
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I. INTRODUCTION 



The recent studies of parity violation in compound nuclear states and the observation 
[ of large longitudinal asymmetries in the transmission of polarized epithermal neutrons 
has lead to a renewed interest in both the parity violating nuclear interaction and the 
relevant structure of the compound nucleus. A considerable number of theoretical papers 
was published on this subject in the last 3-4 years PHi"2f. Most of the theoretical effort 
was directed towards the understanding of the observed enhancements in the asymmetries 
and towards an understanding of sign correlations in these asymmetries detected in the 
studies with 232 Th [§]. The common basis of the majority of these theories was the use of 
a one-body parity non-conserving potential that induces parity mixing in the nucleus. In 
spite of this common basis the various works used different mathematical approaches and 
calculated a number of different terms that contribute to the parity violating asymmetries. 
In the present work, starting from a one-body PV potential, we shall develop a consistent 
framework in which we attempt to calculate or estimate all the contributions emerging from 
such an approach. 

The recent theoretical work as well as the pioneering papers in this field 0|| deal with a 
new type of experiments that were first performed in the early eighties |I[ and then expanded 
in the in the early nineties mainly at the Los Alamos Neutron Scattering Center (LANSCE) 
||§|. In these experiments one observes resonance scattering to p-wave resonances. Using 
polarized neutron beams and unpolarized medium and heavy mass targets one measures the 
longitudinal asymmetries: 

r ^> a+(E)+<j-(E)' { > 

where <r + , a are the total cross sections for neutrons with positive and negative helicities 
respectively. A non zero value for P indicates that parity is not conserved in the process. 
The large enhancements by factors 10 5 predicted by theory were indeed observed in the 
original experiments about ten years ago M as well as in the more recent ones [0,E§. The 
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great advantage of the recent improved Los Alamos experiments is that one is now able to 
measure parity violation of a number of pi/ 2 resonances in the same nucleus. This enabled 
one to apply statistical arguments in order to determine the mean square PV matrix element 
and the average asymmetry P. 

We should stress here that the excitation energies in the neutron plus target nucleus that 
are reached with eV neutrons are of the order of 5 — 10 MeV. In heavy nuclei this is the 
regime of the compound nucleus where the density of states is very large, the spacing between 
for example the J = | + states is of the order of 10 — 20 eV. This means that the wave 
function of each such state contains about 10 6 "principal" components, i.e. components that 
have comparable amplitudes each of the order 10 -3 . The components are defined as many- 
body shell-model configurations resulting from an appropriate single-particle basis. Such 
large complexity of the compound nuclear states of course means that statistical arguments 
should apply. It is therefore quite a surprise that in the case of 232 Th the sign of the 
asymmetry P measured for all the nine states, that show statistically significant results, is 
positive. The average asymmetry for this nucleus was determined || to be about 8%. 

The resonances, both the s and p, are well separated, their width being of the order of 
a very small fraction of an eV, while as mentioned the separation is around 10 eV. Most 
of these resonances were observed in the past in neutron scattering reactions. Higher I 
resonances (d, f etc.) are not seen because of the low penetrabilities of such resonances. 
Among the p-wave resonances that are observed, some have total spins of J = |. These will 
not show up in the asymmetry measurements because parity violation in such resonances 
occurs through the mixing with the d 3 / 2 state and so the low penetrability will hinder the 
observation of any asymmetry signal. 

In the approach taken in this paper, as well as in some other recent theoretical works we 
are able to incorporate aspects of relatively simple nuclear structure and dynamics as well as 
statistical behaviour. In the present work this duality will be quite explicit. In some cases 
we shall average over certain class of compound states while other states will be treated 
individually. 
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II. ASYMMETRY P AND THE T-MATRIX 



The central quantity that is of interest in these parity violating studies is the asymmetry 
given in Eq. (JI|). In all the experiments we discussed the asymmetry is measured on and 
off resonance. The on resonance measurements are for the J = | to which we shall refer 
as the pi/2 resonances. (It should be clear that in compound states the single particle pi/2 
component is very small of the order of 10 -5 admixture). For a.pi/2 resonance the asymmetry 
is: 

P(E) - a+ ^ {E) - a '^ {E) ( 2 ) 

Pl/2 V ' 

where a pi (E) and a pi (E) are the helicity plus and minus on resonance total pi/2 cross 
sections, and a pi (E) is the total helicity independent cross section for such resonance. 
In Appendix [X| we derive the relation between the asymmetry and the T-matrix. Here 
we present this well known || formula relating P(E) to the parity violating part of the 
T-matrix, T PV : 

III. REACTION THEORY 



We shall use the projection theory of Feshbach [13| and its extensions [|14|T5ll to formulate 
the reaction theory for the parity violating experiments in the compound nucleus. Let us 
divide the complete space into the following subspaces relevant to our problem. The {P} 
space (the open channel space) will contain the channels 4> Pl/2 \0) and Sl/2 |O), where |0) is 
the ground state of the target, <fi Pl/2 and <fi Sl , 2 are the continuum single-particle states. The 
corresponding projection operator that projects onto this space will be denoted as P (not to 
be confused with P{E) the asymmetry in Eq. ([!])). The {Q} space will contain all J 71 = | + 
bound states \q) and the projection operator will be denoted as Q. In the {R} space we 
shall keep the J n = | bound states |r) and the corresponding projection operator will be 

5 



denoted by R. The rest of the existing states will not be of interest to us and therefore are 
not included in our considerations. As usual we shall assume that the spaces {P, Q, R} are 
mutually orthogonal, P ■ Q =P ■ R =Q ■ R = and that P 2 = P, Q 2 = Q, R 2 = R. A state 
of total spin J = ~ is then written in terms of the above projection operators as: 

$ = P^ + Q^j + R^j . (4) 

Denoting by H the full nuclear Hamiltonian which contains both parity conserving and 
parity non-conserving parts, we can write the following set of coupled equations: 

(E - H PP )Py = H PR m + Hp Q QV , 

(E - H RR )RV = H RP Pm + HrqQ* , (5) 
(E - H QQ )QV = H Q pPV + H QR RV , 

where the usual notation Hab stands for AHB. We shall now eliminate the explicit reference 
to the part of the wave function by introducing the operator TC, given by: 

H = H + HQ 1 QH . (6) 
& - Uqq 

Then we can write: 

(E - H PP )Py = HprRV , 

(E - H RR )RV = HrpPV . (7) 

The projection of the operator TC in Eq. @ onto the P-space TCpp we shall term as the 
optical model operator. Solving formally one can write for P^\ 

(E - H P p - H PR j? 1 H R p)P^ = . (8) 

Let $( ±} be the solutions of (E - TLpp)^ = , i.e. 

(E - (H + HQ 1 QH) PP )^ = . (9) 

We note that in zero order of the PV interaction the optical model operator contains for 
the p-wave channel only Hpp, while for the s-wave the above operator contains also the more 
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complicated term which involves the |g)-states. We should call the attention to the fact that 
following the experimental approach we attempt to calculate the various observables (cross 
section, asymmetry) for individual p-wave resonances only. The coupled equations (|7|) will 
give rise to rapidly fluctuating scattering amplitudes due to the {Q} space. The averaging 
over this space will be made at a later stage (see Sec. [V]). 



A. Casel. No channel coupling. 

Let us first consider the case when the direct parity violating coupling between the chan- 
nels is neglected. The Lippmann-Schwinger equation for P^f^ (subscript c denotes p or s) 
can be written: 

P*(±) = $(±) 



& — ripp & — rtRR 



(10) 



and more explicitly as: 



P«r(±) = $(±) j 7Y PR 7i R p$ (±) 11) 

& — rtpp & — rtRR — rLRPyti — ripp) ripR 



The T-matrix is: 



± sp — J-sp 



| (€ , \nPR\r)(r\n R p\^ p +> ) 
h'E- (r\H RR \r) - {r\n R pG^n PR \r) 



(12) 



where G P = — — } and T^p is the non-resonant part of the T-matrix due to Hpp. 

From now on we shall consider only the parity- violating part of the T-matrix: 

T S P P V = T^ PV + T^ PV . (13) 



Here Xjr es ) py is the resonant part and T^p PV is due to the parity-violating part of TC 



pp. 



where ($ s ^^p*') is the parity- violating S-matrix. Unless we take channel coupling into 
account, $ c has a definite parity and the parity- violating part of Tjp PV is zero. To shorten 
our notation we write: 

(r\H RR \r) + {r\H RP G { p7ip B \r) = E rT . (15) 

To first order in the parity violating interaction the second term in Eq. fll2|) gives two 
contributions. When H R p is parity non-conserving we get the term: 

t pv _^ (€ ) \V PV \r)(r\H RP \& p +) ) 

sp ' dir v e - E r + |r r ' { ! 

Note that this term although resonant does not involve the mixing with \q) states and 



therefore similarly to the case of isobaric analog resonances |TJ[] we will refer to it as the 
"direct" term. 

When the parity violation appears in RHQ QHP the corresponding part of the 

E - H QQ 

T-matrix, called the compound term, will be: 

t pv _^ (€ ) \Hp Q \q)(q\V PV \r)(r\H RP \<l> ( p +) ) 
sp ' cornp (E — E r + |r r )(£ l — E q ) ' 1 ) 

We should comment here that the various compound states \q) have very small (but nonzero) 
decay widths and therefore the energies E q do contain small imaginary parts preventing the 
expression in Eq. ( p|) and in all forthcoming equations involving E — E q in the denominator 
to became singular. For the sake of simplicity of notation we have not included explicitly 
the above imaginary part. 

B. Case2. Channel coupling. 

We consider now the case when there is direct parity violating coupling between the two 
channels that have opposite parity. We write equation (E — Tlpp)$^ = in the form 



S 



( np + P n + p - P U 



, H PP H PP ^ 



E 



(18) 



or 



(E - if PP - if 



++ 



E-Hoo HqP HpP 



H 



pp 



~ H PQ 

x($ p +$ s ; 



-h; 



^E-Hqq^ 



where $ p and $ s are wave functions of pi/2 and S1/2 channels. The superscripts denote 
parities. Parity- violating terms oiTipp will give the T^p PV part of the T-matrix: 



T (P)PV 
sp 



^^\Hp + Q \q)(q\V PV \K +) ) 



E — E n 



(19) 



To display the channel coupling contribution in the resonant part of the T-matrix we write 
the solution $ for the p— wave channel in first order of the PV-interaction: 

* = % + e^h^ pp + H PQ E -H QQ H « p) * a ' (20) 
where $ c are solutions of the uncoupled equations with definite parity, used in the previous 
section. Now inserting $^ ' for $ s ) in the resonant part of the T-matrix in Eq. (|T2"| ) in the 
case, when TCpp and TLrp are parity-conserving operators, i.e. equal to H PR and iffjp we 
get in the first order of the PV-interaction two terms, resulting from the parity- admixed exit 
channels: 

(C1^^ TT ^/fpH|r)(r|/fpp|< ) ) 



rpPV 

sp,(res,cc) 



= E 



E (+)-H pp - 

E — E r + |T r 



£ 

r,q 



Where: 



E(±) - Hp P 
(subscript cc denotes channel coupling 



(^^Hqp^^V^^^Hpp^^Hpp^) 
{E-E r + \T r ){E-E q ) 



E / dE'\^(E')}^-^(^(E') 



(21) 



(22) 



IV. EXPRESSIONS FOR THE ASYMMETRIES 



To calculate the asymmetry P(E) we make use of Eq. (^) and the equation for the total 
P1/2 cross section given by: 

^4 7 r 2 7 r|(r|^p|< ) )| 2 r r 

where T r is the total width of a pi/2-wave resonance. The statistical weight factor equals 
unity for pu 2 wave |L(J. The on- resonance (i? = T r ) asymmetry resulting from the direct 
term in Eq. (|T6| ) is: 







r 


)(r 


Hrp 




K 


r 




p 


-))| 2 





p^(E r ) = -2 " « — 1 : ' p • (24) 



The on-resonance contribution of the compound part of the T-matrix in Eq. flTT| ) to the 
asymmetry is: 



r J 





H PQ \q)(q\V PV 


r)(r 


Hrp 




\(r 


\Hrp 




-E q ) 





Pcomp\E> 

g IV \i-iRP\v p l\- ■ y-i^r - ^q) 

(25) 

The asymmetry resulting from the channel coupling contribution in the first term of Eq. 

1 

EM - Hp P 



l|) is analogous to the direct term with the matrix element ($g ) \V PV —, . — — Hpp\r) 



replacing (<& s |V py |r) and given by: 



xP.V. / dE 



,(^\E r )\V pv \%{E'))(%{E>)\Hp R \r) 



E r — E' 



(26) 



(The symbol P.V. stands for "principal value"). This contribution was suggested by 
Lewenkopf and Weidenmiiller as a possible source of enhancement in P CC (E). The 



second term of the T-matrix in Eq. fl2l|) gives the following contribution: 
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P cc (E r ) = -2 ' 



x 



n v 



i(ri^pi$r)i 2 

($ri^PQk)(r|^p|<'(^)) 



xP.V. ftfg«W>| . (27) 

Finally we shall write the non- resonant channel coupling contributions to P(E). Because we 
consider here only the pi/2-resonance case we shall evaluate the non- resonant contribution 
at the energy E pa E r meaning that we can neglect the contribution of the tails of s-wave 
resonances. (In the off-resonance region <t° has normally a contribution from the tail of the 
nearest s-resonance, which is usually bigger than the p-wave off-resonance cross section). The 
first term of the non-resonant channel coupling part of the T-matrix Eq. (|i9|) for E = E r 
gives: 

This term is equal to zero, if the the radial parts of cj) p and <j) s functions are purely real. The 
second term in Eq. ( jT9|) yields for E = E r : 

Pcc{Er) = \(r\H RP W;>w 

q E r — E q 

We should stress here that a number of contributions we have written here will be 
numerically small and not have an impact on the final results. The numerical calculations 



and estimates of these various contributions to P(E r ) will be discussed in Sec. [IX. 
One can write all the expressions for the asymmetry in a concise form: 

P(E r ) = -2 Ke{1 ;\ } , (30) 

where 7J is the neutron escape amplitude: 

7r T = (r\H RP \& p +) ) (31) 
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and 7^ py ^ is the parity violating escape amplitude relevant to each of the processes discussed 
above. For example in the direct process, described by Eq. ( |16D the corresponding 7^ py ^ is: 



7<Z ,T = (<'\V r \r) , (32) 



while for the compound process in Eq. (|T 



!r,comp / y jj\ ' V 00 / 



V. DOORWAY STATE APPROACH 

The expressions for the asymmetry P comp (E) in Eq. (|25| ) are very complicated because 
they involve in principle an infinite sum over all compound states This complexity is of 
course inherent in the physics of the compound nucleus, and usually statistical methods are 
adopted in order to calculate the observables P,p7|. However, it was recognized in the past 
that due to the absence of strong multi-particle forces in the nuclear Hamiltonian one can 



often reduce the complexity of the problem by introducing the notion of doorways fL4| , |15 
The idea is to select in the {Q} space states which have direct couplings via a one or two- 
body force either to the channel space or the states in the {R} space. These special selected 
states are the doorways. Formally, we divide the {Q} space into the subspace of doorways 
denoted as {D}, the states belonging to it as \d), and the corresponding projection operator 
as D. The rest of the {Q} space will be denoted as {Q'} and its states as \q'). 

{Q} = {D} + {Q'} (34) 

The states \d) are not eigenstates even of Hqq. The doorway states thus are model states 
that are treated explicitly. The influence of the \q') is taken into account indirectly by 
introducing new parameters into the theory, such as widths, energy shifts, etc. which are 
usually treated phenomenologically. For example in our discussion of Eq. ( p5|) and the PV 
matrix element (r\V PV \q) we replace the \q) states by some doorways \d). These doorways 
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however couple to the \q') states through a strong parity conserving interaction giving rise 
to a width of the doorway. 

At this stage we introduce an average over states \q'). Since we deal in the present 
work with doorways that are of one-particle nature the averaging over \q') states in the 
effective Hamiltonian Eq. (|J) will be the same as the averaging which leads to the appearance 
of a smoothly varying optical model potential. The averaging interval I should satisfy 
D sp > I 3> D q , where D sp and D q are the average spacings between the single-particle states 
and the compound states correspondingly. As a result of this averaging the s-wave channel 
wave functions <3> s will be solutions of a smooth, complex optical model potential Tip p. The 
p-wave channel wave function is not affected by this average and $ p = $ p . 

Let us now concentrate on Eq. (|17]) for the T-matrix and apply to it the doorway state 
approximation. Using now $ s and taking a Lorentzian average of the part containing the 
| q) states we arrive at an expression: 

fPV v (^- ) |gp Q | g >(g|V py |r)(r|gflp|< ) > 

s P ,com P {E _ ^ + t _ v ^ E _ Hqq + i /} • I ) 



We use now the following relation [14.1 



D l -^D = 

E - Hqq + ^1 

E — H DD + |7 — H DQ /[E — Hqiqi + ^I\~ x Hqi d 



Assuming random signs for the matrix elements (q\HQip,\d) one may write [ 14 , 1^ ,Q : 



(36) 



E-Hqq + ^I Y E-E d + \T d 
with 



Ed -iri= m \ d)+i: mmm (38) 

so that Eq. (|1^) becomes: 



2 d ^ 1 ' Y E - E * + ¥ 
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1 sp,comp 2-^i ( jp P , i r \/ r P , irh ' \ ) 



- { E-E r + ^ 7 ){E-E d + {T l d ) 
Using now Eq. (|^) and Eq. (^) we find for the compound contribution to the asymmetry 
at resonance energy E = E r the expression: 

Re{{^;\H PD \d){d\V pv \r){r\H RP \^ ) )){E r -E d ) 



Pcompi^Er) 2 ^ ] 



|2 
4 J 



M^lgp^ldXdl^lrXrlffflpI^))^ (4Q) 



Taking the case when \E r — E d \ we find the expression: 

Pcompi^Er) 





tf Pr ,|d><d|V PV |r><r| 










$; +) )i 2 (^- 







(41) 



Note that this expression can be obtained from Eq. ( p5|) by replacing the states |g) with the 
doorways \d) and the wave function ' with a smooth optical wave function $ s . At this 
stage we have not yet specified the doorways. These will depend on the nature of the H we 
use. 

One can write Eq. © again in a concise form by making use of Eq. ([3(]) and the 
equation: 

where j d (E r ) is the escape amplitude of the doorway: 

yl(E r ) = (&- ) \H PD \d) • (43) 

Eq. ( P7| ) can be treated in the doorway state approximation by replacing the states \q) 
by a doorway. However, the matrix element (d\V PV \<§> p (E')) will be much smaller, than the 
matrix element (d\V PV \r) because the doorway \d) is not related through V FV to <& P (E'). 
We expect therefore that the contribution in Eq. Q27| ) will be a small correction to P comp (E r ) 
in Eq. fl2|). 
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VI. ONE-BODY PARITY VIOLATING POTENTIAL AND THE SPIN-DIPOLE 

DOORWAY 



In the preceeding sections we have discussed in general terms the derivation of the longi- 
tudinal asymmetry without specifying the nature of the parity violating part of the Hamil- 
tonian. Therefore it was not possible to specify the nature of the doorways that might be 
relevant in the calculation of the asymmetry or other quantities related to the experiments. 
We shall now proceed by limiting our treatment to a one-body parity violating potential. 
(In a forthcoming publication we shall address the question of two-body parity violating 
interactions and the question of the relevant doorways in such a case). The one-body poten- 
tial we shall now treat is assumed to originate from a two-body parity violating interaction 



of the type discussed in Ref. |0|. The two-body PV force as discussed in the above ref- 
erence or other work is constructed from meson-exchange models and therefore contains 
a substantial amount of phenomeno logical input. The interactions so constructed do not 
of course incorporate the very short range behaviour of the PV interaction. Convoluting 
such two-body interactions with the nuclear density one obtains a one-body parity violating 
potential. Requiring that in first order such potential will be even under time reversal one 
arrives at a potential that has the form [20[|: 



V PV = elO- 7 ^{/(r i ),^-Pic} , (44) 

where a is the nucleon spin and p = —ih'V the nucleon momentum operator, f(r) is a 
radial function and e is a parameter that designates the strength of the potential. The curly 
brackets denote the anticommutation operation. The weak interaction does not conserve 
isospin therefore the parity violating potential may have an isoscalar and isovector part. 
It was demonstrated |7,^l[| that the doorway that produces a coupling between the state 



r 



) and the \q) states is the spin-dipole giant state |22| built on |r). This doorway is obtained 



by acting with the 

V>; • n ~ Y>[<x, ® Yi^r (45) 

i i 
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operator on |r): 

Ro^Enk^^^rk) , (46) 

where Nq is a normalization constant, Y\ is the vector spherical harmonic of order 1. The 
spin operator is treated as a vector and when coupled to the spherical harmonic Y"i can 
produce an operator that carries total spin J = 0~, 1~, 2~. The state in Eq. fl46|) is of 
isoscalar type and when we multiply the operator in Eq. (^) by the isospin operator t z and 
act on the state |r) we obtain the isovector spin-dipole giant state: 

\SD) 1 = ^-Y, r i[<Ti®Y 1 (r i ))°-T zi \r) , (47) 

iVi 

where N\ is again a normalization constant. We should mention here that in fact it is 
the isovector spin dipole that was studied quite extensively both experimentally p2| , p3| and 



theoretically |p4}| . The information that exists about the isoscalar spin-dipole is rather scarce 
|[25|| . In the work [7[ it was pointed out that the spin-dipole giant resonance exhausts most 
of the strength associated with the a • p operator of the PV potential. (In fact for harmonic 
oscillator wave functions this spin-dipole state exhausts all the strength). The states in Eqs. 
( PS|) and PT| ) are to be used in the expressions in Eq. (pEOp. The dipole giant resonance 
based on the g.s. is located at the position E ~ 78A" 1 / 3 MeV. Experiments |26] show that 



the energy position of a giant dipole resonance built on an excited state is not very different 
from the one built on the g.s. Experimental studies indicate |2^j23|| that the spin-dipole is 
close in energy to the spin independent dipole, therefore: 

E d -E r ~ 78A~ 1/3 MeV . (48) 

The width of the spin-dipole is of the order of a few MeV and thus the condition 
\E d — E r \ > T d is valid and Eq. (f41~D should be used. We are dealing with the limit of distant 
doorways. (In a future publication we shall address the question of nearby doorways). 

We should note the following: in an odd-even nucleus the spin dipole operator when 
acting say on a single-particle pi/2 component of the |r) state will produce a component 
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which will be a s.p. su^ state. This component will be only a part of the \SD) doorway, 
in addition there will be other components which are particle-hole excitations with respect 
to state \r). (When we say "single-particle" we mean a single-particle state with respect 
to the target g.s.). In our present approach as well as the one in Ref. [f7|,|P2| the doorway 
is chosen in such a way as to maximize the matrix element (r|V|rf). We note however that 
in the expression in Eq. fl4l"l) appears in the numerator also the escape amplitude of the 



doorway. The choice of the doorway given in Eqs. ( flq , [47l) does not involve the consideration 
to maximize this quantity. In an approach taken in Ref. || the s\/2 s.p. state is chosen as 
a doorway. This choice does maximize 7^ but does not give consideration to the size of the 
PV matrix element. Both approaches give the same result for P comp (E r ). Since the final 
expression was derived in detail in Ref. [12| we shall present here only the resulting expression 



for the asymmetry P comp (E r ). (For reasons of completeness we repeat the derivation in 
Appendix |B|). The bound s.p. pyz belonging to the {R} space will be denoted as \4>p l/2 )', 
the bound s.p. si/ 2 belonging to the {Q} space will be denoted as \(j>f 1/2 ) and the s.p. escape 
amplitudes will be denoted as: 

7j(^) = « /2 l^p|< +) ); 

1 l{E r ) = & ) \H PD \^J . (49) 

The asymmetry is: 



where Esd,t is an energy, representing the average position of the isoscalar and isovector 
spin-dipole, built on the state |r). We should stress that this equation is derived under the 
assumption, that f(r) in Eq. (|I4]) is a constant. The use of a realistic f(r) will numerically 
changes the result very little. 

The same kind of expression was derived in Ref. || where the doorway as mentioned 
was the state |0^ /2 )- The only difference is that in the denominator of Eq. ( pUp Esd,t is the 
energy position of a collective spin-dipole excitation while in Ref. |§ the energy position of 
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the single-particle state \<j>f 1/2 ) appears. The expression in Eq. (p0[) depends very weakly on 
the particular structure of the state |r). The only reference to \r) is in the denominator, 
but E r — Esd.t has only a very slight dependence on \r). Of course there is a rather strong 
dependence on the energy E r (when E r is small) in the ratio: 

ll(E r ) 1 



4{E r ) VE r 
Note, that the sign of Eq. (|50| ) is fixed for all |r). 



(51) 



VII. PARITY VIOLATING SPREADING WIDTHS 

In the case of model states, thus states that are not eigenstates of the full Hamiltonian 
it is useful to use the notation of a spreading width. This width is a measure of the extent 
to which the model state is not an eigenstate and is spread among the eigenstates of the 
system. If a state is a solution of a Schrodinger equation with a model Hamiltonian having 
some symmetry then one can introduce the coupling via a symmetry breaking potential of 
this state to the rest of space and see how this state is spread out. In this case one can 
speak about a symmetry breaking spreading width. In our case of course it is the reflection 
symmetry which is relevant. We may ask the question what is the parity violating spreading 
width of the negative parity state in the {R} space when coupled to the other two spaces. 
In fact the division into spaces is designed so that the coupling to the {Q} space is only 
through the parity violating potential. The spreading width for any state \r) is given by 

IM!IIMII : 



^-^E^J, (52) 

where I is the averaging interval over the \n) states, such that / ^> D n , where D n is the 
average spacing between states \n). The sum extends in principle over all states, different 
from |r). We shall now consider only the parity- violating part in H (V PV ) and neglect the 
coupling to the {P} space. Eq. ( |52| ) becomes: 



18 



r^-^E#»r. (53) 



q E r — E q + |J 



we 



We may now apply the doorway state approximation to Y\ PV . Following Ref. [f?|,|T^,[T5 
write: 

rr = £ ]{r ™^ d , (54) 
d (E r -E d y + T -f- 

where the sum is over doorways \d) and is the spreading width of the doorway resulting 
from the strong interaction coupling of \d) to the \q') states. The expression for is given 
by 

n=- 2 M E Ji« , (55) 

q> h d ~ E q , + 2 1 

In Fig. H] a schematic representation of the doorway state approach to T d PV is shown. In 
the case of the spin-dipole doorways, we have two, the isoscalar and isovector. Therefore: 

(E - E r f + \ (Et- E r f + 

where Eq, E\, Tq, V[ are the energy centroids and spreading widths of isoscalar and isovector 
spin-dipole giant resonances. In our numerical estimates we shall replace the two doorways 
by one which effectively will represent the two doorways, thus: 

rr = l{rlvPV]s % n , (57, 

(E d - E r ) 2 + V -f 

where Ed, T d is the average energy centroid, and spreading width of the combined distribu- 
tion of isoscalar and isovector spin-dipole strength. The spreading width has the property 
that is weakly dependent on the density of states, consequently also on the excitation en- 
ergy. It can provide therefore a measure of symmetry breaking. As we shall see in Sec. [TX 



the parity violating spreading width was determined in the recent experiments [0,0 and a 
comparison between our estimates and the experimental Tl PV will provide valuable insight 
into the possible mechanisms of parity mixing in the compound nucleus. 
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VIII. ORTHOGONALIZATION 



When we calculate wave functions belonging to the {P} space of open channels and the 
wave functions of the {Q} or {R} spaces, consisting of bound states we find these by solving 
the Schrodinger equation with different Hamiltonians. This leads to the loss of orthogo- 
nality between the above spaces. In order to correct for this one needs to orthogonalize 
the corresponding wave functions explicitly. We do this by orthogonalizing the continuum 
single-particle wave function <fi Pl/2 or 4> Sl/2 to the corresponding single-particle bound state 
wave functions found in the {Q} or {R} spaces. 



The orthogonalization procedure used here is taken from ||15|| , where it was applied to 
the isobaric-analog resonances. Let h be the optical Hamiltonian and |0) the function that 
solves the Schrodinger equation: 

h\4 w ) = E\4 W ) . (58) 

Let us denote by \u) the normalized nuclear single-particle wave function belonging to {R} 
or {Q} space. Then we can write a projected optical Hamiltonian h as: 

h= (1- \u)(u\)h{l- \u)(u\) . (59) 

The projected wave function |0) satisfies the equations: 

h\$ +) ) = E\$ +) ) and (4>\u) = . (60) 

To obtain a solution of these equations we use the following method [ 151 ; after solving 



the homogenous equation (|58D we solve the inhomogeneous equation: 

(E - h)\f) = -\u) , (61) 
then the function is given by: 

\r)=xL w )-\f)S$ff). (62) 
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Here the constant A is determined by requirement of the usual asymptotic normalization of 

In Sec. |X| we shall discuss the numerical results and the effect of the orthogonalization 
on the results. 

IX. CALCULATIONS AND ESTIMATES 

In this section we shall evaluate numerically (or make order of magnitude estimates in 



the case of small terms) the various contributions to the asymmetry P(E r ) derived in Sec. |V 
as well as the PV spreading width discussed in Sec. |VTJ . 

Before we deal with each of the terms separately we present the general ingredients used 
in our calculations. We shall perform all the calculations for the compound nucleus of 233 Th, 
considering it as a typical case. Most of the results that will be calculated for this nucleus 
should apply also to other nuclei in this region of the periodic table, such as U etc. There is 
however one peculiarity to 233 Th, namely that all the nine experimental asymmetries that 
are statistically significant turn out to be positive in this nucleus ||. This kind of sign 
correlation was not found in the other nuclei for which experiments were performed. We 
shall return to this point later. 

A. PV interaction and optical potential 

For the effective one-body PV interaction we used the potential in Eq. QUI) . The radial 
function f(r) has a Woods-Saxon shape with the same parameters as f a (r) used in the 
optical potential Eq. ( |63|) below. 

The continuum wave functions were calculated using the Perey-Buck [B7[ and Madland- 



Young pIpH] optical potentials, and the computer code DWUCK4 [E9f. Both potentials 



have the form: 



Vopt = Vof a {r) + V ls rf s ^ls + iW a b ^ , (63) 
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where /j(r) = [1 + exp((r — i?)/aj)] _1 is the Woods-Saxon function, R = tqA 1 ^. The pa- 
rameters of the Perey-Buck potential are: V = -48.0 MeV, V ts = 14.8 MeV, W = 44 
MeV, ro = ri s = 1.27 fm, a a = ai s = 0.65 fm, a b = 0.47 fm. The parameters of the 
Madland- Young potential are: V = -(50.378 - 27.073^ - 0.345E lab ) MeV, V ls = 24.3 
MeV, W a b = 4(9.265 - 12.666^ - 0.232E lab ) MeV-fm , r = 1.264 fm, r ls = 1.01 fm, 
a a = 0.612 fm, a b = 0.553 + 0.0144:Ei ab fm, ai s = 0.75 fm. Here Ei ab is the neutron en- 
ergy in MeV. The depth parameters of volume, spin-orbit and imaginary surface po- 
tentials for 232 Th and E[ ab = 1 eV are: the Perey-Buck potential, Vq = —48.0 MeV, 
V ls r? s = 23.9 MeVfm 2 , W b = 20.7 MeV-fm; the Madland- Young potential, V = -44.3 
MeV, Vi s rf s = 24.8 MeV-fm 2 , W b = 14.2 MeV-fm. 



In our formalism the optical model operator (Sec. |T|) has a different structure in the 
case of the s and p channels, which is reflected in fact, that we do not average over the 
P1/2, but average over the si/ 2 resonances. Because of that in our calculations of the pi/ 2 
continuum wave functions we do not include the imaginary part of the optical potential. In 
our studies we found, however, that while the asymmetry P(E) is sensitive to the depth of 
the real part of the optical potential, it is not sensitive to the presence or absence of the 
imaginary part of the potential in either the s or p-wave channels. 

The relevant single-particle bound states in the Th region are the 4p!/ 2 , 4si/2, 5si/ 2 - In 
Sec. [VI] they were denoted as \(fip 1/2 ) and \<f>f 1/2 )- The wave functions of these states were 
calculated in a Woods-Saxon potential: 

V sm = V l f s {r) + V^rl- d ^ls , (64) 

r ar 

where a s = 0.60 fm, r s = 1.29 fm, £ = —0.804 and V\ = —69.7 MeV. This set is standard in 



calculations for nonspherical nuclei such as Th ||30|| . This potential produces a deeply bound 
5su2 state, which in nature is in fact unbound. The reason is that when this potential is 
used in deformed nuclei, additional terms are added. These terms are not included in the 
present calculations. To make the state 5sx/2 weakly bound we also used a second set with 
V 1 = -59.5 MeV. 
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In our calculations we checked the sensitivity of the asymmetry P(E) to the depth 
parameters of the optical and shell model potentials to verify the stability of the numerical 
results. 

The strong interaction parts of the Hamiltonians Hpp and Hqp contain both one and 
two-body parts. However as a consequence of the fact, that in the present approach we 
study only the one-body PV potential the contributions of the two-body parts of Hpp and 
Hqp mostly drop out. In some cases the two-body will contribute. This is however beyond 
the scope of the present work and it will be discussed in forthcoming work. 

We use for these two interactions a real Woods-Saxon potential and denote it V s t r - For 
convenience the parameters of these transition potentials are taken from the real part of 
the of the Perey-Buck potential. We should remark that these are not optical potentials 
which are used in the calculation of the channels that are in the {P} space. We should 
also emphasize that in our approach the one-body potentials are not represented by self- 
consistent Hartree-Fock potentials because for example the {R} space contains only bound 
single-particle components while the unbound components are in the {P} space which is 
orthogonal to {R}. This is inherent in the nature of the projection operator approach 

ignis . 

B. Effect of orthogonalization 

When we orthogonalize the continuum wave function (w.f.) to the single-particle bound 
wave function, we can see, that the number of nodes and, hence, the phase at the origin of 
the orthogonalized function is determined by the behaviour of the bound wave function. (Of 
course the asymptotic behaviour of the orthogonalized wave function should not be changed). 
The behaviour of orthogonalized and non-orthogonalized continuum radial wave functions 
is illustrated in Fig. |2|. The continuum s and p w.f., when orthogonalized to the bound 
4px/2 and 5si/ 2 radial wave functions change the phase at the origin, because they acquire 
an additional node. The s-wave continuum w.f., orthogonalized to the 4si/2 radial w.f. is 



23 



just distorted so, that the spacing between nodes is compressed. This explains the changing 
of the phases of the single-particle amplitudes. The strong effect of orthogonalization in the 
present calculations is due to the fact that the continuum wave functions are very close to 



threshold. (This is quite different from the case of analog resonances |Rj| ) . 

Because we are considering here neutrons with energies E n ~ 1 — 500 eV, which are 
much smaller than the depth of the nuclear potential, the behaviour of the continuum wave 
functions inside a nucleus is determined by the potential and is weakly dependent on the 
neutron energy. All the terms, which contribute to P(E) have the same energy dependence 
1 / \[E due to the factor 1 / j p . 

We shall consider as representative the calculation with orthogonalized continuum wave 
functions in the Perey-Buck optical potential and with bound states calculated in the po- 
tential with Vi = -69.7 MeV. 



C. Calculation of P CO mp(E r ) 

We calculate P comp (E r ) in the doorway state approximation. The formula in Eq. (|50"D 
for the on-resonance asymmetry P comp (E r ) involves the matrix element of V PV between s.p. 
bound states and two s.p. escape amplitudes. Our calculation gives for the matrix elements 
{^Pi/i\V PV \nsy2) (in units of e) — iYJ.2 eV and —H7.0 eV for the deeply bound 4s and 
5s states respectively and — il5.8, — zl0.3 eV for the more weakly bound. We used here 
the phase convention of positive radial wave functions near the origin. The corresponding 
matrix elements for oscillator wave functions are — zl6.1 eV and — zl5.2 eV. Using Eq. 
we obtain: 



Pcompi^Pr 





Vstr 


4s V j 


){As 1/2 \V PV \Ap l/2 )(Ap 1/2 
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(65) 



This equation represents the fact, that two doorway states with 5sy 2 an d 4si/ 2 single- 
particle w.f. can contribute. The denominators of these two parts have different signs. In our 
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calculation we used for the energy difference in the denominator the values of ±5 MeV, which 
is supposed to represent the centroid of the isovector and isoscalar spin-dipole doorways. 
The resulting P CO mp(E r ) as a function of E r is shown in Fig. |3|. 



For the orthogonalized continuum waves both terms in Eq. (|65| ) have negative sign and 
the quantity P comp {E T = leV) je — —0.95%. We found in our computations that the presence 
or absence of the imaginary part of the optical potential in the s-wave or in the p-wave chan- 
nel does not change the results noticeably. We also checked the sensitivity of P comp (E r ) to the 
optical and shell model potentials and the orthogonalization of continuum wave functions. 
The asymmetry, calculated with the orthogonalized continuum w.f. is sensitive to the details 
of bound wave functions, thus to the depth of the shell model potential. When V\ is changing 
from —69.7 to —59.5 MeV the asymmetry becomes about 1.5 times smaller for 4s term and 

about 4 times smaller for 5s term and the final result is: P comp (E r = leV)/e = —0.44%. The 

7s 

sensitivity to the depth of the optical potential is less, than 10%. The ratio I — I of single- 

1p 

particle escape amplitudes varies between (0.35 — 0.73) x 10 3 , depending on the number of 

,7s, _. VS , 1o3 



nodes and shell model parameters. The square well estimate is I — I ~ - — = 1 x 10 JT6J. 

7 P kR 



For the non-orthogonalized continuum waves the two terms in Eq. ( |6"5] ) have different 
signs as opposed to the case of orthogonalized waves. The asymmetry, calculated with non- 
orthogonalized waves, is very sensitive to the depth of real part of the optical potential. 
The quantity P comp (E r = leV)/e = 0.06%, when calculated with the Perey-Buck potential 
(V = -48.0 MeV) and 0.34%, when the Madland- Young potential (V = -44.3 MeV) was 



used. The sensitivity to the depth parameter of the shell-model potential is smaller, about 

20%. The ratio 1 — 1= 0.06 - 0.67 x 10 3 . 

1p 

D. Calculation of Pdir(E r ) 

Assuming, that only the one-body part of V PV or Hpr is contributing and writing the 
wave function 



r) = a r \4p 1/2 ) + 
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(66) 



we obtain from Eq. fl24|) that 



P dir (E r ) = -2 





V PV \Ap l/2 )(A Pl/2 


Vstr 
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The resulting asymmetry P<nr{E r ) as a function of the neutron energy is shown in Fig. |j. 
For E r = 1 eV, Pdi r {E r )/e = 0.16%. We see, that the direct term is smaller than the 
compound term by about a factor of 6. The dependence on the depth parameters and 
orthogonalization is of the same type, as found in the case of the compound (doorway) term. 
When V\ = —59.5 MeV, Pdi r {E r = leV)/e = 0.12%. For the non-orthogonalized continuum 
wave function the asymmetry has a negative sign and Pdi r {E r )/e = —0.015% in case of the 
Perey-Buck potential and —0.05% for the Madland- Young potential. 



E. Channel coupling resonance terms 



The term in Eq. (|26"D, as was shown in |]TT|, contains the enhancement, caused by off-shell 
contributions E' 3> E r , because at small energies (j) p (E,r) ~ y/kjl{kr) and the integrand is 
small for E' = E r . The energy dependence of the asymmetry P cc (E r ) is also \ j\[E~ r . For 
E r = 1 eV, P cc (E r )/e = —0.15%, while the unenhanced on-shell result is 1 x 10~ 10 %. This 
asymmetry is not very sensitive to the parameters of the optical potential but shows sensitiv- 
ity to the depth of the shell model potential: P cc (E r )/e = -0.32%, when V x = -59.5 MeV. 
For the non-orthogonalized continuum w.f. the asymmetry P(E r )/e = —0.094% in case of 
the Perey-Buck optical potential and —0.07%, when the Madland- Young potential with a 
constant Vq = —44.3 MeV was used. 



It was suggested in [TT|] to use a surface-peaked approximation, when dealing with inte- 



grals of the kind, found in Eq. ([26|). The assumption is, that the strong and parity- violating 
interactions are surface-peaked and that the equations 

(%{E')\H PR \r) _ (<S> p (E)\H RP \r) 



(j) p {k'R) <p P i kR ) 

and 



(68) 
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^\E r )\V pv \%{E')) _ {^\E r )\V PV \%m 



(j) p (k'R) (j)p(kR) 
hold through the entire virtual energy scale. With the approximation E r = (Ref. UTTp: 

P cc (E r ) = 2Re{(&-\E r )\V PV \%{E r .))) 

r dE' \ME>)\ 2 
Jo E> |0 p (£ r )| 2 ' 

The enhancement here arises due to the off-shell, effect, i.e. when E' 3> 0. 

In our calculations we found, that the main contribution to the integral comes from the 
interval E' = 0.1 — 100 MeV (see Fig. |5](a)) both in the exact calculation and also when the 
above approximation is applied. However, because of the rather complex dependence of the 
continuum wave function (p p (E',r) on energy, the approximations of Eq. (|68|j69|) are good 
only in the energy interval — 1 MeV in the case of orthogonalized and — 0.015 MeV in 
the case of the non-orthogonalized continuum wave functions. In Fig. |5](b) we show the 
ratio of the product of the left hand side of Eqs. (|68| ) and (^) of the product of the right 
hand side of these equations for the case of E = 1 eV as a function of the virtual energy E' . 
We can see that for E' smaller than 1 MeV the ratio is approximately one, and thus the 
surface-peaked approximation works well. For higher E' energies the approximation breaks 
down. Note that the very sharp structure is due to the use of a logarithmic energy scale. 
This conclusion does not change when the derivative of a Woods-Saxon function —f' a {r), 
which is surface-peaked, was used for the strong interaction H RP . 

The value of the asymmetry in the case of the orthogonalized continuum wave functions 
and surface-peaked approximation for E r = 1 eV is P cc (E r )/e = —0.17%, which is close to 
the value obtained in the exact calculation. The agreement might be fortuitous due to 
the oscillations of the exact integrand. This result is not very sensitive to the depth pa- 
rameter of the optical potential. The calculation in the surface-peaked approximation with 
non-orthogonalized continuum wave functions gives 0.61%. As opposed to the previous or- 
thogonalized continuum case, the latter result is very sensitive to the depth of the optical 



potential if V > 46 MeV. In Ref. fy] the value of P cc (E r = leV)/e = 0.25 - 0.75% is found 
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for this asymmetry term. 

The second channel coupling resonance term in Eq. fl2"T|) will be small because of the 
argument given at the end of Sec. |V[ 



F. Channel coupling background terms 



As pointed out in Sec. [TV] the term in Eq. ( Pq) is equal to zero in the limit of real radial 

Im(J) s _ 2 Imtpp _ q 

continuum wave functions cp s and cp v . Because — — — ~ 5 x 10 and — — — ~ 5 x 10 this 

Re<p s Reipp 



term is much smaller, than the one in Eq. (Pq) . The same is true for the term in Eq. (\29\ 
which is much smaller, than the term in Eq. 



We summarize all the results in Table |[ The values P(E r = leV)/e, calculated for 
different terms are shown for the case of orthogonalized continuum wave functions and for 
the two shell-model potentials used in our calculation of the bound states. 



G. Estimate of the parity violating spreading width 

We use Eq. fl57D to estimate the parity violating spreading width in 233 Th. The 
matrix element (r\V PV \d) where \d) is the collective spin-dipole when compared to a 
single-particle matrix elements it is enhanced by a factor y/~N, where N is roughly 
the number of particle-hole states, composing the spin-dipole. In the case of A ~ 230, 
counting the number of particle-hole states, N ~ 100 and we therefore can use the 
estimate (r\V PV \SD) ~ 10 x (<f) p \V PV \4>f ). For e = 1 this gives a matrix element 
(r\V PV \SD) ~ 150 eV. Let us now take E d - E r ~ 5 MeV and ~ 3 MeV assuming that 
this number represents a spreading width common to the isoscalar and isovector spin-dipole. 
A typical spreading width of a giant resonance is of this size. As for the energy difference 
the dipole state is at an excitation energy of about 10 MeV, however we expect the isoscalar 
spin-dipole to be lower and therefore we chose the above value for Ed — E r . We should em- 
phasize however that these numbers are used only in an order of magnitude estimate that 
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we make here. We obtain for T\. PV ~ 2 x 10~ 3 eV, three orders of magnitude larger than 
the experimental value || of T\, PV ~ 6 x 10 -7 eV. In order to obtain agreement one must 
reduce e to a value of about 0.02. Although the above estimate is rough it clearly indicates 
that there is an internal inconsistency in the evaluation of the PV spreading width and aver- 
age asymmetry. The same PV matrix element cannot reproduce simultaneously the average 
asymmetry and the spreading width. The later requires a parity violating potential that 
has a value of e which is "reasonable" compared to the values derived using meson-exchange 



models [IS 



X. DISCUSSION OF RESULTS AND CONCLUSIONS 

The calculations of the average asymmetry P(E r ) in the framework of the one-body 
approximation and for realistic e = 0.1 — 0.3 give numbers which are of the order of 0.1% or 
less. We have considered several contributions, some new, as for example Pd ir (E r ). We have 
evaluated the channel coupling term, suggested in Ref. [[□]], exactly, avoiding the surface- 



peaked approximation. This term turned out to be of the order of 0.02%, if the above 
value of e is used. Even if all these contributions add coherently the asymmetry found, 
would not exceed 0.2%. The average asymmetry found in the experiments for Th and other 
nuclei is of the order of several percent, thus at least one order of magnitude smaller than 
the present theoretical results. Moreover, when using our doorway model to estimate the 
spreading width we find that the matrix element required to reproduce the experimental 
PV spreading width is of the order of a fraction of an eV which corresponds to a value 
of e = 0.02. Saying it in other words, if we use the value of e that fits the asymmetry 
P(E r ) we shall overestimate the spreading width by about two orders of magnitude. We feel 
that we have exhausted most of the contributions to P(E r ) in the framework of one-body 
theories. The above difficulty to reproduce the correct order of magnitude of the average 
asymmetry with reasonable values for e and the internal contradiction within the model 
between the estimate of the spreading width and the calculation of P(E r ) leads us to the 
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conclusion that the one-particle model is not sufficient in trying to account for all the aspects 
of parity violation effects in the compound nucleus. One should attempt to go beyond the 
one-particle approximation and try to introduce the effects of the residual two-body force. 



In that case questions of intermediate structure |H,5H in the nucleus become relevant. In 
a future publication we shall address this question. 

It was emphasized in the literature f||,§,|Il],|l^] that the measured asymmetries for the 



individual resonances exhibit sign correlations in 233 Th. When averaged over all the observed 
resonances |r) the averaged P(E r ) turns out to be around 8%. We have seen that most of 
the terms we considered here, the doorway state approximation to the compound term, the 
direct term, as well as the channel coupling term, have signs independent on the particular 
compound resonance |r). The problem of course is that the size of these terms for the 
realistic values of the PV coupling constant is, as said above, at least an order of magnitude 
smaller than the value found in the 233 Th experiment ||. 

In summary, in this work using a unified reaction theory approach we have derived in a 
systematic way all the terms that contribute to the non-fluctuating part of the longitudinal 
asymmetry P{E). Off-resonance and on- resonance terms were identified. All these terms 
were treated numerically either by direct calculation or in the case of small terms, using 
estimates. The projection operator reaction theory we used is flexible enough so that by 
choosing appropriately the subspaces of the full space we were able to reduce the number 
of terms in the description of the asymmetry. The compound term emerged as the most 
important one. This term was calculated using the doorway state approximation with the 
doorway being the giant spin-dipole resonance. 
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APPENDIX A: 
T-MATRIX AND ASYMMETRY P(E) 



The asymmetry P(E) is defined as: 

P(E) - 2°IJE) • (A1) 
where + and — denote helicities u of projectile neutrons. In this Appendix we derive the 
relation between the T-matrix and P(E), given in Eq. (|3|). Writing the optical theorem for 
helicity v : 

%Imf vu {k^k) = Ev><7w (A2) 

we get: 

° + - v~ = £(Jm/ ++ (0) - 7m/__(0)) . (A3) 
The partial wave expansion for helicity amplitudes f v > v (0) is p^j : 

M0) = E j (2j + l)4'u(0)ft'u, (A4) 

where fi, u are partial amplitudes. In the case when j = \ only and 6 = 0: / ++ (0) = 2f^^ 2 , 
/__(0) = 2/i= 1/2 and 

a + - <T = f (Jm/|; 1/2 - Jm/i= 1/2 ) . (A5) 
We go now to the Isj representation, assuming rotational invariance: 



\Ejmjv) = £ J$k(lOl"\jv)\Els = \ ]3 m 3 ) . (A6) 
i 

The expression for helicity partial amplitudes in terms of partial amplitudes in the Isj 
representation fl x x = f 3 v l is: 

1 2 ' 2 



= E vt ; (^o^bv><ro^bv>/^ tl . (A?) 

in 2 ' 2 

Here 1,1' = j ± g and 
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/i + = K/U^+4^- 2 -^-i>' < A8) 

Because of time-reversal invariance /// 1 = fl v . Then for the case of j = | : 

a + - a' = 8 f(Imf{= + 1/2 - Imf=y 2 ) 

= -If l m f izi l=0 . (AS) 

Expressing the T-matrix elements in terms of partial amplitudes T 3 VI = —~fp ; we obtain: 

1 } ~ ^fjE) " ~W (E) • (A10) 

APPENDIX B: 
DERIVATION OF EQUATION (M) FOR P(E r ) 



In Eqs. ([|I|,[|^,[|7]) we specified the isoscalar and isovector doorways \d) as: 

l^)o = ^E(^|r), (Bl) 
7V o i=l 

In what follows we shall consider one doorway state \d) for each |r) with the normalization 
constant A r . For the sake of simplicity we shall also take the one-body parity-violating 

A 

potential in Eq. ([44]) in its schematic form V PV = g^(<xp)j with g being a constant. The 

i=i 

product of matrix elements of V PV and H^p operators in the numerator of Eqs. fl39,41) is 



now: 



(& s -'\H PD \d){d\V pv \r) 

A A 

= 4<€~^M><r|E(^£(^k> • (B2) 
iV i=i i=i 

The matrix elements of the operator crp can be approximately replaced (for example, using 

the basis of oscillator wave functions) by the matrix elements of iumcrr (see for example 
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Ref. (T^j). (In a recent microscopic calculation |33|] the distribution of erp and err strengths 



was found to be very close to each other). Therefore: 

A A A A 

(r\ XX CTr )i J2( a P)i\ r ) = iujm ( r \ H( crr )i H( fTr )il r ) • ( B3 ) 

i=l i=l i=l i=l 

The last matrix element in this equation is the total spin-dipole strength and therefore equals 
to N 2 . Eq. ([|D is now: 

(&; ] \H PD \d)(d\V PV \r) = gioomN^^Hp^d) . (B4) 

Let us write the state |r) in the form: 

|r> = ar|< /a > + |/> (B5) 

where |/) stands for more complicated components. The wave function of the spin-dipole 
one can write as: 

\SD) = ^« /2 kr-|< i/2 )a r |0f i/2 ) + \f) , (B6) 

where |/') again are complicated states. Because we take into account only the one-body 
part of Hop only the single-particle component in the spin-dipole doorway will contribute. 
Thus: 

(^{Hpn^V^lr) 

= ^maXC'l^l^)^!^!^) • (B7) 
One can write for harmonic oscillator wave functions: 

gi^Jar^J = « /2 |^|< 1/2 ) (B8) 

The last matrix element in Eq. (|B7| ) is just the single-particle escape amplitude "fl(E s = E T ) 
in Eq. (^9|). Using one-body potentials, only the single-particle part a>r\4>p 1/2 ) will contribute 



1/2 

P 



to the matrix element (r\HpR\<& ) in Eq. fl^T|). The coefficients a r will drop out and we 
finally obtain: 

( \4{E r )\*.{E r -E d ) 
where jl(E r ) is the singe-particle escape amplitude for the pi/ 2 resonance, given in Eq. (|49|). 
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FIGURES 

FIG. 1. A schematic representation of the doorway state approximation in the calculation of 
the parity violating spreading width (see text). 

FIG. 2. (a) The unprojected and projected p 1 / 2 continuum wave functions for E n = 1 eV. 
(b) The unprojected and projected S1/2 continuum wave functions for E n = 1 eV. The wave func- 
tions are given in arbitrary units. 

FIG. 3. The calculated compound term of the asymmetry P CO mp using the doorway state ap- 
proximation. The results are given as a function of the neutron energy E n . The two curves 
correspond to calculations in which different shell model potentials were used. 

FIG. 4. The calculated direct term of the asymmetry P calculated as a function of the neutron 
energy E n . The two curves correspond to calculations in which different shell model potentials 
were used. 



FIG. 5. (a) The integral in Eq. ( |70|) calculated as a function of the upper limit virtual energy 
E' (arbitrary units), (b) A comparison between the exact and surface-peaked approximation. The 
ratio of the product of the left hand side of Eqs. (pq) and (69) of the product of the right hand 
sides of these equations is shown for the case of E = 1 eV as a function of the virtual energy E'. 
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TABLES 

TABLE I. The values P(E r = leV)/e, for different terms for the case of orthogonalized con- 
tinuum wave functions calculated for two shell model potentials. Potential I has Vq = —69.7 MeV, 

potential II has V = -59.5 MeV. 

P{E r = leV)/e in % 

term potential I potential II 

compound (doorway) -0.95 -0.44 

direct 0.16 0.12 

channel coupling resonance -0.32 -0.15 

J2Pi -111 -0.47 
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